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I. INTRODUCTION 


Theories featuring gauge bosons, fermions, and scalars constitute the backbone of the standard model of particle 
interactions. This is to date one of the most successful models of nature. The recent discovery of asymptotically safe 
quantum field theories in four space-time dimensions [1] , including their quantum-corrected potentials [2] , widens the 
horizon of fundamental theories that can be used beyond the traditional asymptotically free paradigm [31 [3] . The 
novelty resides in the occurrence of an exact interacting ultraviolet (UV) fixed point rather than a UV non-interacting 
fixed point, as it is the case for asymptotically free theories. 

For the class of theories we will be investigating here a crucial property was unveiled in Ref. [1]: the Yukawa 
interactions, mediated by the scalars, compensate for the loss of asymptotic freedom due to the large number of 
gauged fermion flavours and therefore cure the subsequent growth of the gauge coupling. The further interplay of the 
gauge, Yukawa, and scalar interactions ensures that all couplings reach a stable interacting UV hxed point allowing 
for a complete asymptotic safety scenario in all couplings [1] . This is different from the complete asymptotic freedom 
scenario liHZj where all couplings vanish in the UV, see Refs. mm for recent studies. 

The phase diagram including the scaling exponents of the theory was determined to the maximum known order in 
perturbation theory mm- It was also shown in Ref. [3] that the scalar potential is stable at the classical and quantum 
level. Therefore these theories hold a special status: at arbitrarily short scales and without assuming additional 
symmetries they are fundamental according to Wilson’s definition. 

Having at our disposal new classes of fundamental theories one can use them to construct new dark-matter paradigms 
[To] or even support cosmic inflation HU. It is therefore timely as well as theoretically and phenomenologically 
relevant to investigate, in a controllable manner, the thermodynamics of four-dimensional completely asymptotically 
safe theorie^ 

We organise this paper as follows. The theory and its salient zero-temperature properties are reviewed in Sec. [l^ 
This is followed by the determination of the asymptotically safe pressure to the leading order (TO), next-to-leading 
(NLO), and next-to-next-to-leading order (NNLO) in Sec. Ill The entropy density of the system is determined and its 
properties are discussed in Sec. EYi The thermal degrees of freedom count for asymptotically safe theories is introduced 
and discussed in Sec. |V[ We offer our conclusions in Sec. VI where we also briefly discuss the impact of introducing 
a quark chemical potential. In the appendix we report the beta functions of the theory and further details of the 
computations of the pressure at the respective orders in perturbation theory. 


II. ZERO-TEMPERATURE PHYSICS 

Here we briefly review the salient aspects of the gauge-Yukawa system introduced in Ref. (U such as the phase 
diagram of the theory and the expressions for the UV-safe trajectories away from the UV-stable hxed point. We will 
also provide the expressions for the running of the couplings along the globally dehned UV-IR connecting line known 


^ The thermodynamics of asymptotically free theories featuring perturbatively controllable and interacting IR fixed points has been 
investigated in Ref. m to the maximum known order in perturbation theory. 
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as separatrix, or line of physics. Further quantities required for the subsequent thermodynamical analysis will also be 
reported. 

The asymptotically safe theory suggested in Ref. [T] contains — 1 non-Abelian gauge fields, N^Nj massless 
Dirac fermions, ■(/;, and N'j massless complex-valued scalars, H. The theory has a local SU{Nc) gauge symmetry and 
a global U{Nf)L x U{Nf)ii chiral symmetry at the classical level. Because of the Adler-Bell-Jackiw axial anomaly 
the quantum global symmetry is SU{Nf)L x SU{Nf)ii x U{l)v- The left- and right-handed fermions live in the 
fundamental {Nf,0) and {0,Nf) representations of this symmetry group, 

1pL - = UL-tpL , = UriPh , ( 1 ) 

while the scalars live in the adjoint {Nf,N^) representation, 

H ^ H' = UlHUI. (2) 

It is convenient to decompose the complex-valued (TV/ x TV/) matrix H in terms of the generators Ta in the fundamental 
representation of U (TV/), a = 0,..., TV| — 1, 

H = {Sa + iPa)Ta , (3) 

where Sa are TV| scalar and Pa are TV| pseudoscalar fields. The Lagrangian reads 

^ = -\ V' + Tr {d^H^ d^H) +y{^LHP^R + - uTt {H^Pf - v [Tr {H^H)]" , (4) 

with the covariant derivative 


P>ti=d^ + . 


( 5 ) 


Here, g is the coupling constant of the non-Abelian gauge sector and T^, i = 1,...,TV^ — 1, are the generators of 
SU{Nc) in the fundamental representation. 

At the classical level the theory counts four marginal couplings, the gauge coupling g, the Yukawa coupling y, the 
quartic scalar coupling u, and the double-trace scalar coupling v, which we write as: 


g^Na y^Na uNf v Nj 

( 4^)2 ’ ( 4^)2 ’ ( 4^)2 ’ ( 4^)2 ■ 


( 6 ) 


The appropriate powers of TVc and TV/ in the normalization of the couplings allow to take the Veneziano limit of 
the theory. Following Ref. [1] we will also use the short-hand notation /3i = dtai, with i = {g,y,h,v), for the beta 
functions of the respective couplings (§. It is convenient to introduce the continuous real parameter 


Ni_n 

Na 2 


( 7 ) 


in the Veneziano limit of large TV/ and TVc, with the ratio Nf/Nc fixed. The relevant beta functions of the theory 
have been obtained in Ref. m in dimensional regularisation, using the results of Refs. |14H16j . and are summarised 
in App. I^in the Veneziano limit. 

For e < 0 the theory is asymptotically free in the gauge sector while for e > 0 it becomes a non-Abelian QED-like 
theory because asymptotic freedom is lost. It was shown in Ref. [T] that in this latter case the theory exhibits an 
interacting UV fixed point in all four couplings. This fixed point is controllable in perturbation theory, provided 
0 < e ^ 1. The existence of such an interacting UV fixed point ensures that the theory is a fundamental one, i.e., it is 
valid at arbitrarily short and large distances. Furthermore, the scalar interactions are free from the triviality problem 
because of the presence of an interacting UV fixed point. Therefore, the elementary scalars are part of a Wilsonian 
fundamental theory. 

After a lengthy study of the zeros of the theory [1] , using the gauge-matter system ( |Al| - ( |A4| ), and including also 
the investigation of the stability of the associated classical and quantum scalar potentia [5] one arrives at the only 
mathematically and physically acceptable fixed point accessible in perturbation theory: 


r.* - 26 , I 23(75245-13068^23) ,2 , 

57 370386 ^ ) 

* _ 4 : ( 43549 2300U^ \ ,2 , 3\ 

^2/ “ 19 '"^1^20577 6859 J ^ J 

< = ^e + Oie^), 

q:*i = (2v^ - ^20 -f 6v^) e -f 0{e^) . 
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FIG. 1. Review of the phase diagram of Ref. mm- The gauge-Yukawa subsector of couplings (ag,ay) is shown to leading 
order in the left-hand panel and the gauge-Yukawa-scalar subsector {ag,ay,ah) at next-to-next-to-leading order accuracy in 
the right-hand panel for e = 0.05. Also shown are the UV and IR fixed points (dots), the UV-safe trajectories (thick red line). 
A few trajectories are highlighted as thin magenta lines, and a few generic trajectories are shown as thin gray lines. Arrows 
point towards the IR. Further details can be found in Ref. mm- 


Here we give the analytic expression of the fixed point in an expansion in the small e parameter. 

The phase diagram of the theory was established in Ref. [T] at next-to-leading order accuracy and extended to 
the next-to-next-to-leading order in Ref. [2] where the effects from the running scalar couplings were considered. In 
order to keep the paper self-contained we summarise in Fig. [r the phase diagram of the theory shown in Ref. [5]. 
In the left-hand panel we show the RG trajectories for the Jag, ay) couplings, while in the right-hand panel the 
three-dimensional RG flow is illustrated, that includes also the coupling ah- The two plots include the UV and IR 
fixed points. We have also indicated in the left panel the relevant and irrelevant directions dictated by the signs of 
the scaling exponents. 

The IR fixed point is non-interacting and it is therefore located at the origin of coupling space. The thick red line 
connects the IR and UV fixed point and therefore is the UV-complete trajectory that we term the line of physics. 
The line of physics is also known as separatrix since it separates different regions of the theory in RG space. Along 
the line of physics the theory is non-interacting in the deep IR. The separatrix continues beyond the UV fixed point 
towards large couplings in the IR, leading to a strongly coupled theory presumably breaking conformality and chiral 
symmetry in the IR. The region of the RG phase diagram emanating from the UV-stable fixed point and leading to 
stable trajectories is known as the UV critical surface. Here this critical surface is one-dimensional [T] and it has a 
dynamical nature. 

As mentioned above, the separatrix connects the UV fixed point with the Gaussian one and it agrees with the UV 
critical surface near the fixed point [1] . Although one can always determine numerically the globally defined separatrix, 
it is illuminating, in view of their use in the thermodynamical analysis, to consider an analytical approximation that 
is accurate in the limit of vanishing e. This leads to the following relations among the couplings along the separatrix 

[ 2 ]: 


ay — ^2 , 

~ ^ ag , 

ay = ^ - 2 \/^^ ag . 


(9) 


The one-dimensional nature of the line of physics is encoded in the fact that it is sufhcient to know the running of 
the gauge coupling in order to determine the running of all the other couplings. The precise analytic running of the 
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gauge coupling was determined in Ref. [5] and reads 


?(m) = 


where W{fj,) = is the Lambert function satisfying the relation z = W expfL with 


W= -^-1 

Otn 


and 


2 = 


^ - 1 exp ^ - 1 


( 10 ) 


( 11 ) 


Here is the value of the gauge coupling at the scale /xq, with /x//io ranging between 0 and oo and the gauge coupling 
ranging between 0 < ag < Og- 

Inserting Eq. (10) into Eq. (|^ yields an analytic description of the RG evolution of all couplings along the line of 
physics. This constitutes the zero-temperature information we need to establish the thermodynamical properties of 
the theory. 

At asymptotically high energies W (/i) vanishes while it grows towards the infrared. It is convenient to fix a° via 


Og = a*/(l -|- k) with k G K+, which in practice amounts to fixing the arbitrary renormalization reference scale /xq 
along the RG flow. As pointed out in Ref. [5] the value k = 1/2, i.e., = 2a*/3, corresponds to an exact critical 

transition scale /xq = Ac above which the physics is dominated by the interacting UV fixed point and below which 
it is governed by the Gaussian IR fixed point. The interacting nature of the UV fixed point is expressed by the fact 
that it is approached as a power law in the renormalization scale 


104^2 

ag{^^):z,a; + iaya;)(^£^ , ( 12 ) 

where fio = /xo(l -|- ©(e)). We have used Eq. and that, in the deep-UV limit, the Lambert function approaches 
zero as 

104^2 

lim W{fi) (X (—] . (13) 

Ai/Mo-foo \k-oJ 


There are several nice and distinctive features of the analytically controllable and completely asymptotically safe 
dynamics presented above. In particular, it constitutes an ideal laboratory to investigate thermodynamical properties 
of the theory that for certain aspects resembles Af = 4 theory. One of the similarities is the fact that along the line 
of physics all couplings are related. This is also a basic feature of A/” = 4 theory due, however, to the high degree 
of space-time supersymmetry. In the non-supersymmetric case the relations among the couplings are dynamical in 
nature being dictated by the dimension of the critical surface. In Eig. [^we show in the left panel the beta function 
of the gauge coupling along the line of physics linking the Gaussian fixed point with the interacting UV fixed point. 
In the right panel the running of all couplings along the line of physics is shown using Eqs. ([^ and ( [To| ). From the 




aia Log(f;//Jo) 


FIG. 2. Left panel: Gauge beta function along the line of physics displaying the characteristic asymptotically safe behaviour. 
Right panel: From top to bottom the running of the gauge, Yukawa, single trace, single trace plus double trace, and only single 
trace coupling is shown along the line of physics. We have chosen fe = 1/2 and e = 0.1. 
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figure the completely asymptotically safe nature of the theory is evident. For the sake of completeness we mention 
that, at fixed Nc and large TV/, it has been argued m that an asymptotically safe theory can emerge without scalars 
and to leading order in ^/Nj. Further physical properties of this intriguing possibility were investigated in Ref. [T]. 


III. ASYMPTOTICALLY SAFE THERMODYNAMICS 

We will now study thermodynamical quantities of the theory to leading (LO), next-to-leading (NLO), and next-to- 
next-to-leading order (NNLO) in the couplings. The details of the computation of the thermodynamic pressure for the 
gauge-Yukawa theory, that are valid for any number of colours and flavours and applicable as well to the entire phase 
diagram of the theory, are provided in App. We work here in the Veneziano limit and along the line of physics (|^. 

From Fig. one can immediately see that there are three relevant energy regions with distinct dynamics: the one 
dominated by the Gaussian IR hxed point, i.e., /i ^ Ac, the one dominated by the interacting UV hxed point ^ Ac, 
and the cross-over energy region for which /i ~ Ac. By identifying, for example, the renormalization scale with the 
temperature, at zero chemical potential, we can test the thermodynamic properties of the asymptotically safe plasma 
along the entire line of physics. 


A. Hot Asymptotically Safe Pressure to Leading Order 


Near the Gaussian IR fixed point the theory is non-interacting, and the ideal-gas limit applies. This constitutes the 
LO contribution along the entire line of physics. Specialising the results of App. |^for the pressure to the Veneziano 
limit (Vc, Nf ^ 1) at zero chemical potential, and normalising it to the one of the gluons (in the same limit) we have: 


Po 

Po,g 




7^ 

4 Ac ■ 


(14) 


Here we notice that the ratio depends at most quadratically on Nf/Nc- We are, however, considering an expansion 
in e = Nf/Nc — 11/2 and therefore re-express the result in terms of e: 


Po 

Po,g 


327 

” 8 “ 



(15) 


We note that already to this order the thermodynamical expression depends at most on and that furthermore the 
e expansion allows for a new handle on the thermodynamical expansion, which is absent for a generic gauge-Yukawa 
theory. 


B. Hot Asymptotically Safe Pressure to Next-To-Leading Order 


The previous LO expression for the pressure is exact at the IR fixed point because of its non-interacting nature. Since 
we choose to identify the renormalization scale with the temperature, this then occurs for very small temperatures, 
i.e., n = T Ac- However, when the temperature rises, the plasma starts feeling the various interactions. To NLO 
the pressure as function of the temperature reads: 


Po +2 

Po,g 


327 

“8” 



ag + (ay + 2afi) 



25 

T 



(16) 


Besides the trivial scaling which cancels between numerator and denominator, there is an additional temperature 
dependence due to the running of the couplings which are evaluated at the temperature T. In a conformal field theory, 
however, we can only consider ratios of scales or, equivalently, the couplings are measured in units of a reference value. 
Along the line of physics the natural choice for the reference scale is Ac, corresponding to a value of the gauge coupling, 
which is 2/3 of its fixed-point value. This is the scale above which the physics is dominated by the UV hxed point, 
while below it is governed by the Gaussian IR one. This allows us to immediately determine the two limiting values 
of the pressure obtained for T <C Ac and for T ^ Ac. 

For T Ac the physics is dominated, as already mentioned earlier, by the non-interacting hxed point and therefore 
the NLO pressure coincides with its ideal-gas expression. 


Po +2 

P0,g 


327 


—e £2 = 40.875 -f 12.75 e -f , 


r< Ac. 


(17) 
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However, for T 3> Ac, i.e., near or at the UV fixed point we can use the fixed-point values for the couplings which 
yields: 


Po+2 

Po,g 


= 40.875 - 84.6877e 0{e^) , T > Ac 


(18) 


At zero temperature and chemical potential we have solved the theory to the maximum known order in perturbation 
theory that abides the Weyl consistency conditions. This implies that we know the gauge and Yukawa couplings to 
the second order in e and the scalar couplings to the leading order in e. This limits the expansion of the pressure in 
powers of e to the next-to-leading order in e. Interestingly, we observe a net drop of the pressure when normalized to 
the ideal-gas limit, valid in the deep IR, due to the interacting nature of the asymptotically safe plasma. 

We can also determine the pressure along the entire line of physics by using Eqs. and ( |l0| ), where Eq. ( [T^ 
reduces to 


= 40.875 -b 12.75e + - (213.613 -b 69.6094e -b 5.75995e^)ag(r) , 

Po,g 


for any T . 


We use Eq. (101 and replace the renormalisation scale fi with T and the reference scale with Ac and write: 

^ 1 -b W{T) ■ 


(19) 


( 20 ) 


From the knowledge of the functional dependence on T we deduce that the coefficient of the pressure decreases 
monotonically when the temperature increases. The theory assumes ideal-gas behaviour only in the deep IR. This is 
different from the case of asymptotically free field theories where the ideal-gas limit is approached in the deep UV. 
From Fig. we observe a decrease in the value of the pressure around Ac, now normalised to the ideal-gas limit, when 



Log(T/po) 

FIG. 3. Pressure normalised to the leading-order ideal-gas value, up to NLO corrections, as function of the temperature. We 
have chosen k = 1/2, meaning that fio = Ac. From bottom to top e assumes the values 0.08 (red), 0.07 (blue), and 0.05 (black). 


increasing the temperature. By increasing e the change in the pressure is more pronounced. 


C. Hot Asymptotically Safe Pressure to Next-To-Next-To-Leading Order 

At this order we observe the emergence of non-analytic contributions in the couplings. More specifically, the leading 
contributions will start at 0{g^) and u^/^). These come from the plasmon-ring diagrams CHI, cf. Fig. and 

the detailed computation of their contribution to the pressure is given in App. The respective contribution to the 
pressure in the Veneziano limit reads: 

-= (15 -b 2e)5 -|—— {‘2oih -b ciy -b oiy)^ (11 -b 2e)^ = 150\/5otg -|--b cXy -b cuy)^ -b 0{e^) . (21) 

Po,g v3 v3 v3 
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Here we used the fact that all couplings are already of order e, see Eq. (|^. Along the line of physics this reduces 
to: 


— = 704.061 a„(T) 5+0(65). (22) 

Po,g 


Besides the fact that the contribution is non analytic in e we learn that it starts at and that it is posi¬ 

tive, differently from the NLO contribution which is negative and starts at 0(e). Near the UV fixed point, i.e., at 
temperatures T ^ Ac, it assumes the value: 


— = 216.899 e5 -h ), T > Ac . (23) 

Po.g 

The full pressure to at nonzero temperature, zero chemical potential, in the Veneziano limit, and along the 

line of physics reads: 


PO+2+3 

Po,g 


= 40.875 -f 12.75e - 213.613ag(T) -h 704.061 ag{T)i + 0{e^) , 


for any T . 


(24) 


The pressure normalised to the non-interacting limit is shown in the two panels of Fig. [^for different values of e. For 
the plot on the left the values for e are, from bottom to top, 0.05 (black), 0.03 (blue), and 0.01 (red), while the values 
for the solid curves on the right are 0.08 (red), 0.07 (blue), and again, for reference, 0.05 (black). For the dashed 
curves we retain some of the higher-order corrections in e in order to gain an estimate on their order of magnitude. 
These come by retaining all the powers in e from Eqs. (19) and (21). Overall the analysis shows that we have good 
control of the perturbative expansion up to NNLO terms. Because the NNLO corrections to the pressure are positive. 




FIG. 4. Pressure up to NNLO, normalised to its LO value, as function of the temperature. We have chosen fc = 1/2 meaning 
that /io = Ac. Left panel: The values for e are, from bottom to top, 0.05 (black), 0.03 (blue), and 0.01 (red). Right panel: The 
e values for the solid curves are 0.08 (red), 0.07 (blue), and 0.05 (black). We kept some of the higher-order corrections in t for 
the respective dashed curves. 


the total value of the pressure increases with respect to its NLO value, but not with respect to its LO value, when 
approaching the UV fixed point. Furthermore, the normalised pressure starts developing a minimum near Ac when e 
increases above the value 0.05. 

The decrease of the pressure normalised to the one of the ideal gas is not guaranteed to be monotonous because 
of the effects of the NNLO corrections. However, within the realm of perturbation theory this quantity globally 
decreases, i.e.. 


Ap 


norm 


P0+2+3{T <C Ac) 
Po 


P0+2+3{T > Ac) 

Po 


= 2.384e- 5.306e5 > 0 , 


(25) 


provided that e < 0.202. This is guaranteed by the fact that the radius of convergence of the expansion is, at zero 
temperature, e < 0.11 [T], and it is even smaller at nonzero temperature. 












IV. ASYMPTOTICALLY SAFE ENTROPY 


Another important quantity to determine is the entropy density of the system, which is related to the pressure via 


s = 


dp 
dT ■ 


(26) 


Given that in the present system, and along the line of physics, the pressure can be written as a function of only one 
coupling we have: 

p = /(ag(T))^T4 . (27) 

In the non-interacting gas f{ag(T)) is the number of boson degrees of freedom plus 7/4 times the number of Weyl 
fermions. The entropy density normalised to the one of an ideal gas of gluons reads: 

s 1 

^ = 2(iV2 - 1) 


/ + 


/3(«g) df 

4 da„ 


P , /3(as) d{p/po,g) 


Po,, 


+ 


dan 


P 

Po,g 


1 djp/po^g) 
4 dlnT 


(28) 


with /3{ag) = dag/dlnp,, where p, = T, is the gauge beta function along the line of physics. We have used the fact 
that / = 2(V/ — l)p/po^g. Because the beta functions vanish at a fixed point we have that at the IR and UV fixed 
points the normalised entropy density agrees with the normalised fixed-point pressure and therefore: 


Sir 

S0,g 

SUV 

^0,g 


lim - 

T/An^O Po.g 


fiR 7 Nf , 

If = 1+^ + 4^= 40.875+ 12.75e + e^ 


lim — = ^ = 40.875 - 84.6877e -b 216 . 8996^/2 2 ^ 

T/A,^ooPo.g ^''"2 ‘rAr 2 ^ ^ 


2V2 2V2 


Away from the fixed points the normalised entropy density and pressure differ by the quantity 


(29) 

(30) 


^_ 7 ^ ^ I3{ag) d{p/po^g) 

^0,g PO^g 4 dag 


(31) 


which is directly proportional to the beta function of the theory. This behaviour is different, for example, from the 
case of A/" = 4 theory where the beta function vanishes identically. 

Using the NNLO expression for the pressure from Eq. (24) we deduce 


50+2+3 _ PO+2+3 
50,g P0,g 


-53.4033 


dagjT) 

dlnT 


1 -4.94395 ag(T) 2 


(32) 


This contribution being directly proportional to the gauge beta function along the line of physics from Fig.j^it is clear 
that it is suppressed compared to the term directly proportional to the normalised pressure and therefore decreases 
from the IR to the UV. In Fig. we show the difference in entropy density and pressure both normalised with respect 
to the ideal-gas limit, rather than the ideal gas of gluons. The difference is simply an overall numerical factor. 

We learn that the entropy density normalised to the ideal-gas limit decreases overall from the IR to the UV with the 
dominant piece, in perturbation theory, given by the normalised pressure. However it does not decrease monotonically. 


V. ASYMPTOTICALLY SAFE THERMAL DEGREES OF FREEDOM 

The free energy density J^{T) = —p{T) can be used to count the physical degrees of freedom of the theory at 
different energy scales. The temperature probes the relevant degrees of freedom by exciting them. Therefore the 
function 


f{T) 


F{T) 90 p{T) 90 
r4 )/2 - ’ 


(33) 


is a possible candidate to count these degrees of freedom. Alternatively one can use the coefficient of the opportunely 
normalised entropy density, i.e., / -I- d//4dlnT. The two definitions coincide at fixed points. We can now count, for 
the first time, the thermal degrees of freedom along the entire line of physics of a completely asymptotically safe field 
theory. Up to an overall normalisation the result is the one presented in Fig. It shows that the thermal degrees of 
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FIG. 5. We show the difference between the entropy density normalised to its leading-order ideal-gas value minus the similarly 
normalised pressure, in units of 10“®, up to NNLO corrections as function of the temperature. We have chosen fc = 1/2 meaning 
that fio = Ac and from bottom to top e assumes the values 0.07 (magenta), 0.05 (black), and 0.03 (blue). 


freedom decrease from the infrared to the ultraviolet, albeit not monotonically. In the deep infrared, i.e., for the cold 
field theory, we have of course the ideal gas result 


fjR = lim /(T) = 2(iV2 - 1) + 2N] + , 

T—>0 2 


(34) 


that in the Veneziano limit reads 


2iV2 Ny A ■ 


By construction this function coincides with Eq. (14) and provides the overall normalisation. 
Interestingly we discover 


(35) 


fiR > fuv , 


for the classes of asymptotically safe theories investigated here with 


fuv = lim /(T) . 

T—yoo 


(36) 

(37) 


This demonstrates that the inequality fm < fuv that has been conjectured to be valid for asymptotically free field 
theories [19] does not apply to asymptotically safe field theorie^ The fact that this function decreases in the present 
case is due to the fact that the theory becomes interacting in the UV while it is free from interactions in the deep 
infrared. On the other hand the a-theorem is satished as shown in Ref. m- 


VI. CONCLUSIONS 

We have computed relevant thermodynamic properties of non-supersymmetric four-dimensional completely asymp¬ 
totically safe field theories [Tj up to NNLO. Because of the completely asymptotically safe nature of the theories that 
have been investigated here the coupling constants freeze at a constant and calculable value in the UV. Furthermore, 
their value can be made arbitrarily small in the Veneziano limit because of the existence of a continuous control 
parameter. This has justified a perturbative determination of the vacuum and in-medium properties of the theories 


2 The inequality has not been proven even for asymptotically free theories, but it was shown to be consistent with known results and 
then used to derive constraints for several strongly coupled, vector-like gauge theories. The correct counting of the infrared degrees of 
freedom, with respect to the inequality, for the important case of an SU{2) = Sp{2) gauge theory with fermions in the fundamental 
representation was first performed correctly in Ref. m- The conjecture has been used also for chiral gauge theories |21lI22| . 
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investigated here. In this work we have determined the pressure and entropy density of these theories to next-to- 
next-to-leading order. We find that because of the nature of the interactions, both the pressure and the entropy 
density, normalised to their respective non-interacting ideal-gas values decrease when going from the infrared to the 
ultraviolet. 

After this initial investigation several novel avenues can be explored such as the response of completely asymptoti¬ 
cally safe theories to the introduction of different kinds of chemical potentials. Here, we would just like to remark on 
one interesting physical phenomenon. For the sake of simplicity, we assume that there is a single chemical potential 
fj,q associated to fermion number conservation. Similarly to QCD, at nonzero Hq and sufSciently low temperature, 
the theory features a colour-superconducting phase, because of attractive one-gauge field exchange interactions near 
the Fermi surface [23]. Because of the Pauli principle, fermion Cooper pairs must form in channels which are totally 
antisymmetric in colour-flavour-spin space. To be definite, let us focus on Cooper pairs in the antisymmetric spin- 
zero channel. Such pairs must then be either completely antisymmetric or completely symmetric in both colour and 
flavour. In QCD, the one-gluon exchange interaction is attractive in the antisymmetric colour-antitriplet channel, 
which requires also an antisymmetric wavefunction in flavour space. In principle, this is different in the gauge-Yukawa 
theories studied here: the scalar fields can mediate attractive interactions also for symmetric representations in colour 
space, which in turn demands a symmetric flavour wavefunction. However, whether this actually happens requires 
a more quantitative study, since a repulsive interaction in the symmetric colour channel may destroy the pairing. 
Nevertheless, the phenomenon of colour superconductivity (in an antisymmetric colour channel) in these theories is 
robust: any attractive interaction, no matter how small, will destabilize the Fermi surface and lead to the formation 
of Cooper pairs. Therefore, Cooper pairs will form except right at the Gaussian IR fixed point, i.e., at T = /i = 0. 
In the perturbative regime, i.e., for e ^ 1, the system is a BCS superconductor, with a gap which is exponentially 
small in the coupling. For e <C 1, i.e., Ab ~ 2A^f/ll, a chiral-density wave phase is not expected to occur; this would 
require W > lOOONf Plj. 
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Appendix A: Beta functions 


In the large-A limit, the perturbative renormalisation group equations for the couplings § have been obtained in 
Ref. [T2| in dimensional regularisation, also using the results of Refs. [TH - nT>] . In terms of Eq. Q they are given by 
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(Al) 

(A2) 

(A3) 

(A4) 


for I3g,l3y,l3h, and (3^ up to (3, 2,1, l)-loop order, respectively. In the terminology of Ref. [Tj we refer to this as the 
next-to-next-to-leading order (NNLO) approximation. The NLO approximation corresponds to the approximation 
in which the (2,1,0, 0)-loop terms for I3g,(3y,(3h, and f3y are retained. As discussed in Ref. [1], this ordering of 
perturbation theory is also favoured by the Weyl consistency conditions [131 ESI [26] . 


Appendix B: Explicit computation of the Thermodynamical Properties 

In this appendix we determine the thermodynamical quantities of the theory to leading (LO) and next-to-leading 
order (NLO) in the couplings. For a perturbative calculation of the pressure we need to identify the interaction 
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vertices in Eq. Q. The gauge sector will be analogous to QCD, so we shall simply use the results from Ref. |18j . 
For the Yukawa interaction and the self-interaction of the scalar, however, we will be more explicit. We therefore 
decompose the Yukawa interaction term ~ y in Eq. ffl with the help of Eq. and the definition of the projectors 
onto right- and left-handed chirality, Pb.,l = (1 ± 75)/2, 

i’R Sa + i'4’Ta^5'>P Pa ■ (Bl) 

In order to compute the self-interaction terms ~ rt, n of the scalar field H, we utilize the decomposition (§ , the 
orthogonality relation 


Tr (T,Tb) = ^dab, 

and the (anti-)commutation relations for the generators of U{Nf), 

{Ta, Tb} = dabc Tc , 
[Ta,Tb\ = ifabcTc , 


where dabc {fabc) are the totally (anti-)symmetric structure constants of U{N). We obtain 
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(B5) 


1. Pressure to LO 


To LO, the pressure of the theory Q is that of an ultrarelativistic ideal gas of — 1 gauge fields, NcNf Dirac 
fermions, and ‘^N'j scalars. At temperature T, and if we assume a common chemical potential for all fermions 
(associated to net-fermion number conservation), we have 


where m 


Po{T,Hq) =po,g{T) +poj{T,pq) +po^h{T) , 
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PojiT, Pq) 
Po,h{T) 
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(B6) 


(B7) 


2. Pressure to NLO 


The NLO contribution to the pressure, 

P2{T,Pq) =p 2 ^g{T) P P2,gf{T, Pq) P P2,Hf{T, Pq) +P2,h(7') , 


(B8) 


has a diagrammatic representation in terms of the two-loop diagrams shown in Fig. The first term in Eq. (B8| is 
the contribution from the self-interaction of the gauge fields and gauge fields with ghosts, cf. Figs. I^a-c). It reads 


P2,g{T) 


-g^Na{Nl-l) 


144 ■ 


(B9) 
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FIG. 6. Two-loop contributions to the pressure. Diagrams (a-d) are the same as in QCD; gauge fields are represented by wavy 
lines, ghosts by dotted lines, and fermions by solid lines. Diagrams (e,f) are the two-loop contributions in a Yukawa theory 
where fermions interact with scalar (dashed line) and pseudoscalar fields (dash-dotted line). Diagrams (g-i) are the two-loop 
contributions arising from the self-interaction of the scalar field H, decomposed in terms of scalars and pseudoscalars. The 
combinatorial factors Ci and C 2 are implicitly computed in the text, but will not be needed explicitly. 


The second term in Eg. (|B8[ ) is the contribution from the fermion loop. Fig. j^d), where the fermion interacts with a 
gauge field. This reads |18j 





(BIO) 


The third term in Eg. ( |B8[ ) is the contribution from the fermion loop, Figs.[^e,f), where the fermion interacts either 
with a scalar or a pseudoscalar field. It turns out that both contributions are identical. The calculation proceeds 
analogous to the one for the diagramj^d). The result is, up to a prefactor, identical to Eg. (BIO), 


P2,HfiT,n,) = -y^N]N, + ^ 


Mg 

576 ' 327r2 ' 


(Bll) 


The last term in Eg. (B8) receives contributions from the vertices ^ u and ^ u in Eg. (|^ , 

P 2 ,h{T) = P 2 ,u{T) +P2,viT) ■ 


We first compute the latter. With the help of Eg. (B4) we write 


- V [Tr - ^ (SaSaSbSb + 2SaSaPbPb + PaPaPbPb) 


(B12) 


(B13) 


In order to produce a double-bubble diagram of the type shown in Fig. |^g), either (1) we can tie a leg ^ Sa together 
with the other leg ~ Sa (then we must tie Sb together with Sb), or (2) we can tie a leg ^ Sa together with one of 
the two legs ^ Sb (then the other leg ~ Sa must be tied together with the other leg ^ Sb)- Case (1) corresponds 
to a Hartree-type routing of internal indices and produces an overall factor of Nj, because there are N'j scalar fields 
running in each of the two loops. Case (2) corresponds to a Fock-type routing of internal indices and produces an 
overall factor of N'j, because all indices are tied together in a way that there is effectively only one loop in the index 
a. Overall, we obtain a factor of + 2). The same can be repeated for the pseudoscalar contribution, Fig. ^i), 

with the same result, as nothing distinguishes the two types of fields in the absence of chiral symmetry breaking, xhe 
remaining diagram is the one with one scalar and one pseudoscalar loop, cf. Fig. ih). Here, the factor is (with the 
factor 2 from Eg. (B13l) simply 2N'j. Each double-bubble diagram is proportional to the sguare of a tadpole which, 
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for massless particles, has the value T^/12 [T5]. Altogether we obtain 

P^AT) = -vN%N] + 1)^^. 


(B14) 


We now compute P 2 ui'^ )- We first consider the contribution to the diagrams in Figs, ^g) and (i), corresponding to 
the first line in Eq. (B5). The three different ways to tie legs together to form a double-bubble diagram can be written 
in terms of a combination of Kronecker-deltas, so that the prefactor becomes 

{dabndcdn ^acn dhdn dadndhcn) {^ah^cd ^ac^bd ^ad^bc) — g aan dbbn + “^dabndabn) — ^ + 1) > (B15) 

where we have used daan = V^N^^Sno (a sum over a is implied) and dabndabn = + !)■ (This can be proven 

using the relationship dijkdijk = {N'j — — 4:)/Nf, i,j,k = 1,... ,iV| — 1, for the symmetric structure constants 

of SU{Nf).) Now we consider the contribution to the diagram in Fig.[^h), arising from the second line in Eq. (B5|. 
Here, there is only one way to tie the legs together. 


- {dabndcdn T facnfbdn T fadnfbcn) dabbed — . i,daandccn + 2facnfacn)= l,Nfi2Nj-l) , 


(B16) 


where we have used the relationship fabnfabn = fijkfijk = Nf{N‘j — 1). Putting all this together (remembering that 
Eq. (B15) is multiplied by a factor of 2, for scalar and pseudoscalar contributions), we obtain 


rp4 rp4 

P 2 AT) = -^Nf{2Nj + 1 + 2iV| - 1)— = -2uNf 


144 


144 


(B17) 


Adding the contributions (B9), (BIO), ( Bll[ ), and ( |B12 ) (which is a sum of Eqs. (B14| and ( B17| )), the complete NLO 
contribution to the pressure is 


rTi4 ( iA 

P2 (T, p,) = - [g^ (7V2 - l)iV, + vN] [N] + 1) + 2u fV|] — - [g^ {N^ - l)fV/ + v^N]N,] 
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Using the properly normalised large Nc and Nf couplings and taking the Veneziano limit we arrive at 
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Trading Nf/Nc for e we have: 
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We observe the explicit dependence on the couplings of the theory. 


(B20) 


Appendix C: Pressure to NNLO 

To NNLO, the pressure at nonzero temperature and chemical potential receives contributions from the so-called 
plasmon ring diagrams |18j . cf. Fig. The leading contribution of these plasmon ring diagrams are ~ 0{g^) and 
^ 0{u^d‘^, u^/^). These are calculated in the following. There is a plasmon ring for the gauge field and one for the 
scalar field H, 


P3(T,Pq) = P3^g(T, fig) + P3 ^h(T, Pg) . 

The contribution from the gauge field, Fig. [^a), is the same as in QCD [TB], thus 
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FIG. 7. Plasmon-ring contributions to the pressure. The diagrams (a) are the same as in QCD, while the diagrams (b) are 
the plasmon-ring contribution from the scalar field H (shown only for a scalar Sa, the one for a pseudoscalar Pa has the same 
value). 


with the electric screening mass of the gauge field [18] 


Wg(T,/i,) =5^ 




In analogy to Eq. (C21, the plasmon-ring contribution from the scalar field H reads 


P‘i,H{T,Hq) = 2N‘j 




Un 


where 


= E(0,0) 


(C3) 


(C4) 


(C5) 


is the screening mass of the scalar field; E(0, 0) is the zero-Matsubara frequency, zero-momentum limit of the cor¬ 
responding one-loop self-energy. The prefactor 2iV| takes into account that we have 2iV| (pseudo-)scalar degrees 
of freedom. Thus, we only need to compute E(0, 0) for one of these fields, say the scalar field Sq. In general, the 
one-loop self-energy of the scalar field can be computed by functional differentiation of the two-loop contribution to 
the pressure with respect to the scalar propagator S m, 


E-2 


Sp2 

6S 


(C6) 


This corresponds to amputating a scalar propagator in the diagrams shown in Fig. |6| Obviously, only the diagrams 
(e), (g), and (h) can contribute to E. The contribution from the fermion loop in Fig. ^e) is 


E/(Q) = 2/2t^ 


d^fc 




(27r 


Tr[g{K)TogiK -Q)To] 

npik - ^jLq) - npjp - Pq) npjk -I- Hg) - npjp + Pq) 
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(C7) 


Here, p = k — q, k = k/k, and np{k =F pq) = i jg ^tie Fermi-Dirac distribution for (anti-)particles. 

After renormalization of the vacuum contribution and taking the limit go = 0, (f —>■ 0, we obtain 


S/(0,0) = 
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The contributions from the double-bubble diagrams, Figs, [^g) and (h), can be decomposed into two parts, one 
proportional to the vertex v and one to the vertex u. For the first one, we have from the Hartree-type routing of 
internal indices in Fig. |^g) a factor 2 x Nj (a factor 2 because one can open either one of the two tadpoles), and 
from the Fock-type routing a factor 2x2 = 4 (one factor of 2 because the Fock-type diagram appears twice relative 
to the Hartree-type one and one factor of 2 because one can open either one of the two tadpoles). Finally, we have a 
factor 2x1 = 2 from the diagram in Fig. [^h) (one factor of 2 because this type of diagram appears twice relative to 
Hartree-type one, and a factor of 1 because functional differentiation only opens the scalar tadpole). The remaining 
tadpole is the same, since there is nothing that distinguishes scalar from pseudoscalar fields. After renormalization 
of the tadpole (which then is equal to T^/12) we obtain 




(C9) 


Finally, we compute the contribution from the diagrams in Figs. (g) and (h) proportional to the vertex u. Here we 
go back to the left-hand side of Eq. (B15). Amputating a scalar propagator (corresponding to the zeroth scalar field 


^o) means that any two of the indices a, b, c, d may take the value 0. For symmetry reasons we may restrict ourselves 
to the first factor Sabred in the second set of parentheses and multiply the result by a factor of 3. Now we have either 
a = 6 = 0orc = d = 0. This gives 
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We also need to consider Eq. (B16). Here, only the indices a, b can take the value 0, and thus we have 

^ {dabndcdn 4” facnfbdn 4- fadnfbcn') daO^boSed — ^ doQndccn ~ ^ ■ 
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Putting everything together, we have 


Su = 2u Nf 


T2 

12 


(C12) 


Adding Eqs. (C8), (C9), and (C12|, we finally get 
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